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ON ■ Abstract 

In this paper we provide a computational approach to the shape of curves which are 
u: rational in polar coordinates, i.e. which are defined by means of a parametrization 

U\ \ (r(t),0(t)) where both r(t),9(t) are rational functions. Our study includes theoret- 

ical aspects on the shape of these curves, and algorithmic results which eventually 
lead to an algorithm for plotting the "interesting parts" of the curve, i.e. the parts 



showing the main geometrical features. 



1 Introduction 



Plotting and correct visualization of algebraic curves, both in the case when 
they are implicitly defined by means of a polynomial f(x, y) = 0, or by a 
parametrization (p(t) = (x(t),y(t)) with x(t),y(t) being rational, have re- 
ceived a great deal of attention in the literature on scientific computation 
(see for example pQ, [3], [6], [7], [8], [9], [II]). With this paper, we want to 
initiate a similar study for curves which are written in polar coordinates. Such 
curves may appear in Engineering and Physics, in particular in Mechanics 
(very specially in Celestial Mechanics), and also in Cartography. In this sense, 
here we address those curves which are rational when considered in polar 
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form, i.e. curves defined by means of a parametrization (r(t), 0(t)) where both 
r(t),9(t) are rational functions. Arquimedes' spiral, Cote's spiral under cer- 
tain conditions, Fermat's spiral or Lituus' spiral (some of these can be found 
in www.mathematische-basteleien.de/spiral.htm), for example, belong to this 
class; also, in Cartography these curves may appear when using for example 
conic, pseudo-conic or polyconic projections ([3], [T3]). 

It is quite natural to start wondering if this kind of curves can be algebraic, 
when considered in cartesian coordinates. However, one may prove (see The- 
orem [1] in Subsection 12. ip that, with the exceptions of lines and circles, this 
cannot happen. As a consequence they can exhibit properties that algebraic 
curves cannot have. Here we analyze some of these which are relevant from 
the point of view of plotting, in particular the existence or not of infinitely 
many self-intersections (see Subsection 12. 2p . the appearance of limit points, 
limit circles and spiral branches, all of them introduced in Subsection I2.3[ and 
the relationship between both phenomena. 

The possibility of detecting in advance the above phenomena is very useful to 
improve the plotting of these curves, which can be really "devilish" . To give an 
idea, one may use the computer algebra system Maple, where the instruction 
polarplot is available (we have also tried other packages like Maxima, SAGE 
and Mathematica; however, either no similar instruction was available, or the 
corresponding command behaved in a very similar way to polarplot). This in- 
struction allows to draw curves defined in polar coordinates either by means of 
an equation r = f(6), or by a parametrization (r(t), 0(t)), and has a nice per- 
formance for simple curves, but may not be enough for illustrating the behavior 
of a more complicated curve. As an example, one may consider the following 



ize these curves, one obtains the outputs in Fig. 1: here the curves (1), (2), (3) 
are displayed from left to right; in (1) and (2) we have asked Maple to plot the 
curve for t G (— oo, oo), and in the case of (3) we have chosen t e (—2.01, 2.01) 
(because r, 9 are both non-bounded for t — > ±oo). However, in the three cases 
it is clear that the output does not really help to understand the behavior of 
the curve. 

By using our ideas, one can obtain an algorithm which provides, for a given 
curve of the considered kind, information and several plottings corresponding 
to the more interesting parts of it. We have implemented this algorithm in 
Maple 15. Our algorithm analyzes the curve, detects its main features, and 
uses the command polarplot for plotting the curve over several intervals 
(generated by our algorithm), each one showing certain features of the curve; 
together with the provided information, these plottings help to clarify the 
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Fig. 1. Some complicated curves plotted with polarplot 



behavior of the curve. The outputs of our algorithm for curves in Fig. 1 can 
be checked in Subsection 3.2 and Subsection 3.4. 

The paper is structured as follows: in Section 2, the reader will find the the- 
oretical results on the shape of these curves. In Section 3, we provide some 
details on the algorithm, together with several examples of outputs. Finally, 
in Section 4 we present some conclusions and suggest future lines of research. 



2 Geometrical Properties of Curves Rational in Polar Coordinates 

2. 1 Preliminaries and First Properties 

Along the paper, given a point ? 6 R 2 we denote its polar coordinates by (r, 9), 
and its cartesian coordinates by (x, y). Notice that the polar coordinates of P 
are not unique, since (r,6), (r, 9 + 2kir) with k G Z, or (— r, 9+ (2£+ l)ir) with 
£ G Z, define the same point. 

We will analyze the geometry of a planar curve C which is rational when 
parametrized in polar coordinates. In other words, there exist four real poly- 
nomials A(t), B(t), C(t), D(t) not all of them constant, with gcd(A,B) = 
gcd(C, D) = 1, such that 



parametrizes C in polar coordinates. We will also consider the curve C* = (f(C), 
i.e. the rational planar curve parametrized by ip(t) over the (r, #)-plane; for 
example, if ip(t) = (t,t) (which parametrizes an Arquimedes spiral, then C* is 
a line over the (r, 9) plane. Notice that our goal is not to describe C*, but still 




this curve is useful at several steps for proving certain facts on C. 

Furthermore we assume that <p(t), as a parametrization of C*, is proper, i.e. 
that it is injective for almost all values of t (equivalently, that almost all points 
of C* are generated by just one t- value). Observe that if cp(t) is not proper 
then it can always be properly reparametrized by applying the algorithm in 
Chapter 6 of [12]. Also, we will say that a point Pq G C* is reached by the 
parametrization if there exists some io £ C such that <p(to) = Pq. It can be 
proven (see Proposition 4.2 in [2]) that the only point of C* that may not be 
reached by the parametrization is 

Poo = \im t ^ ±00 <f(t) 

whenever it exists. Furthermore, if <p(t) is proper then the real points of C* 
that are reached by complex, but not real, values of the parameter are, ex- 
cept perhaps for P^, isolated points of C* (see also Proposition 4.2 in [2]). 
Therefore, the real part of C* can be expressed as 

<p(M) U Poo U {Finitely many isolated points.} 

In the sequel we will discard the isolated points of C* and we will focus on 
the remaining part of C*. Since polar coordinates provide a natural mapping 
II (not 1:1) from C* onto C, we will say that a point of C is reached by (p(t) 
if it has the form II(P), where P G C* has been reached by ip(t). In this 
language, notice that we are identifying C = n(C*) where C* is the real part 
of C*, discarding isolated singularities. 

Now if C is either a circle centered at the origin or a line passing through the 
origin, then it is algebraic in algebraic and polar form at the same time (in 
fact, C* is represented by an equation of the type r — r = or 9 — 9 = 0). 
The next theorem proves that circles and lines are the only cases when this 
phenomenon happens. We acknowledge here the help of Fernando San Segundo 
for proving the theorem. Here we use the expression real curve to mean a curve 
with infinitely many real points. 

Theorem 1 Let C be a real algebraic curve. Then it is rational in polar coor- 
dinates if and only if it is either a real circle centered at the origin or a real 
line passing through the origin. 

Proof. The implication (<=) is clear; so, let us focus on (=>■). For this purpose, 
let f(r, 9) G R[r, 9] be the implicit equation of C*. Now since C* is rational, 
then it is irreducible, and so also f(r,9). If / depends only on r, taking into 
account that / is irreducible and that C* is a real curve (because r(t),9(t) 
are real functions), / must be linear. So, C can be either a circle or a point 
(the origin, if fir, 9) = r); but in this last case C is not a real curve, and 
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hence it must be a circle. If / depends only on 9, then for similar reasons 
C must be a line. Thus, assume that / depends on both r, 0, and let us see 
that in this case C cannot be algebraic. Indeed, if C is algebraic, let g(x,y) G 
R[x, y] be its implicit equation and let h(r,9) = g(rcos(0), rsin(0)). Notice 
that every zero of f(r, 9) is also a zero of h(r, 9). In this situation, consider the 
resultant M(9,u,v) = Res r (f(r,9),h(r,u,v)), where h(r,u,v) is the function 
obtained when substituting formally sin(6 l ) = u and cos(0) = v in h(r,9). 
Since C* is irreducible and / depends on both r, 9, does not have any factor 
only depending on r, and therefore f(r, 9) and h(r, u, v) cannot have any 
factor in common; hence, M(9, u, v) cannot be identically 0. Observe also that 
M(9, u, v) must depend explicitly on 9, i.e. M ^ M(u, v). Now by Lemma 4.3.1 
in [H] we have that the resultant Res r (/(r, 9), h(r, u, v)) specializes properly 
when u = sin(6 l ) and v = cos(#), i.e. 

M(0, sin(0), cos(0)) = Res r (/(r, 0), h(r, sin(0), cos(0))). 



Since / depends on both r, we can find an open interval J C M such that 
for every O G I, the equation /(r, O ) = and therefore also h(r,9 ) = 0, 
have at least one real solution. Hence, by well-known properties of resul- 
tants, for all 0o G / it holds that = 0o is a zero of M(0, sin(0), cos(0)), 
i.e. M(0, sin(0), cos(0)) vanishes over /. Hence, by the Identity Theorem (see 
page 81 in [10]), and taking into account that M is analytic, it holds that 
M(0, sin(0), cos(0)) is identically 0. Since M is a polynomial in 0, sin(0), 
cos(0) where is present (i.e. M ^ M(sin(0), cos(0))), this implies that 9, 
sin(0), cos(0), are algebraically dependent. However, this cannot happen be- 
cause sin(0), cos(0) are trascendental functions. □ 

Hereafter we exclude the cases when either r(t) or 9{t) are constant. 



2.2 Self-intersections 



A point generated by t G K corresponds to a self-intersection of C if it exists 
s G M, t s, and k G Z, such that (£, s, k) is solution of some of the following 
two systems: 




r(t) = r(s) 

9{t) = 9(s) + 2k7r 




r(t) = —r(s) 

9(t) = 9(s) + (2k + 1)tt 
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Moreover, when P^ = (r^^oo) exists, we also have to consider the self- 
intersections involving P^, i.e. the solutions of 

. "^oo = r(s) Too = —r(s) 

6 00 = 6{s) + 2kir 0^ = 6{s) + (2k + 1)tt 



So, in order to study the self-intersections of C we have to study the above 
systems. Additionally, whenever the equation r[t) = has more than one 
solution, or one solution and = 0, the origin is also a self- intersection. 

In the algebraic case, Bezout's theorem forces every algebraic curve to have 
finitely many self-intersections. However, in our case this does not necessarily 
hold. Since r(t) is rational, it cannot happen that a same point of C is crossed 
by infinitely many branches of C. But there can be infinitely many different 
self-intersections, which can be detected by analyzing and (*)2- We start 
with (*)i. For this purpose, we write 

r(t)-r(s) = ^ t)B{s) - A{s)m 



8(t) - 9{s) - 2kTT 



B(t)B(s) 

C(t)D(s) - C(s)D(t) - 2k7rD(t)D(s) 
D(t)D(s) 



We denote the numerator of r(t) — r(s) by a(t, s), and the numerator of 9(t) — 
e(s)-2kn by /3(t,s,k). 

Lemma 2 There do not exist a(t, s) ^ 1 and ko £ Z (in fact, ko £ R) such that 
a(t, s) simultaneously divides gcd(a(t, s),(3(t, s, ko)) and B(t)-B(s)-D(t)-D(s). 

Proof. Assume that we have an irreducible polynomial a(t, s) and ko £ Z such 
that a(t,s) divides both gcd(a(t, s), (3(t, s, k )) and B(t) ■ B(s) ■ D(t) ■ D(s). We 
need to show that a(t,s) = 1. Suppose that a(t,s) divides B(t) (resp. B(s)). 
Since it also divides gcd(a(t, s), j3(t, s, k )), it divides a(t, s), and therefore 
it also divides A{t) (resp. A(s)). On the other hand, gcd(A(t), B(t)) = 1 by 
hypothesis. That shows that a(t, s) = 1. If a(t, s) divides D(t) (resp. D(s)) we 
argue similarly with f3(t, s,k). □ 

Then we are ready to proceed with the following theorem, that will have 
important consequences on the study of 

Theorem 3 For any k ^ 0, k £ Z, (in fact, k £ M.) the system has 
finitely many solutions. 
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Proof. In order to prove the statement, we need to show that for k G Z, k ^ 0, 
gcd(a(t, s), j3(t, s, k)) = 1. For this purpose, assume by contradiction that 
there exists some k G Z, k ^ 0, such that H(t,s) = gcd(a(t, s), f3(t, s, k)) ^ 
1. Then H(t,s) defines an algebraic curve H over C 2 ; furthermore, since by 
hypothesis k ^ 0, H(t, s) cannot be t — s. So, there are infinitely many points 
(to, So) G Ji with to ^ Sq. Now by Lemma [2] only finitely many of them 
fulfill B{t) ■ B(s) ■ D(t) ■ D(s) = 0. So, the functions r(t), r(s), 0(t) and 
9(s) are well-defined at almost all the points (t , s ) G H, and r(£ ) = r (so), 
@(to) = @{ s o) + 2A;7r. Since for a fixed to there can only be finitely many values 
of s such that r(t ) = r(s) (because r(t) is rational) we deduce then that there 
are infinitely many points (r ,9 ) G C* such that (ro,0o + 2kn) G C* too. In 
other words, the curves defined over the (r, #)-plane by f(r, 9) and f(r, 6+2kir) 
have infinitely many points in common, and since / is irreducible and both 
have the same degree, then they must define C*. Now let (r , 8q) G C*; then 
( r o> $o + 2/c7r) is a zero of /(r, 6 1 + 2kn), and since /(r, 9 + 2fc7r) also defines C* 
then it is also a point of C*. Following the same reasoning, we conclude that 
(?"0) + 4/c7r) also belongs to C*, and in fact that (r , 6*o + 2nkn) G for all 
nGN. Since ^ we have that these are all different points of C*. Hence, we 
have that C* intersects the line r = r at infinitely many points. But this is 
impossible because C* is algebraic. □ 

We can obtain similar results for (*)2- For this purpose, we denote the numer- 
ator of r(t) + r(s) by fi(t, s), and the numerator of 6{t) — 9(s) — (2k + l)n by 
u(t, s, k). Then the following lemma, analogous to Lemma El holds. 

Lemma 4 There does not exist b(t, s) ^ 1 and k G Z (in fact, ko G M.) such 
that b(t, s) simultaneously divides gcd(/x(t, s), u(t, s, ko)) and B(t) ■ B(s) ■ D(t) ■ 
D(s). 

The following theorem, very similar to Theorem EJ holds. 

Theorem 5 For any k G Z (in fact, ko G R), it holds that gcd(/i(t, s), z/(t, s, /c)) = 
1. 

Proof. Arguing by contradiction as in the proof of Theorem [31 we conclude 
that f(r,9) and /(— r, 9 + 2kix) define the same curve (namely, C*). So, let 
(ro, #o) G C* where r — r does not divide f(r, 9). Then (— r , 9q + (2k + l)ir) G 
C*, and for the same reason (ro, 9q + (4/c + 2)ir) G too. Proceeding this way 
we get that all the points of the form (ro,2n(k + l)7r), with n G N, belong 
to C*. For any value k E Z all these points are different; so, we get that the 
intersection of C* with the line r = r consists of infinitely many different 
points. But this cannot happen because C* is algebraic. □ 

Hence, by Theorem [3] and Theorem [51 we obtain the following result on the 
existence of infinitely many self-intersections of C. 
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Corollary 6 C has infinitely many self-intersections if and only if(*)i or (*)2 
have solutions for infinitely many values k G Z. 

Based on Corollary [61 we have the following result which provides a sufficient 
condition for C to have finitely many self-intersections. 

Theorem 7 If 9(t) is bounded, then there are at most finitely many self- 
intersections. 

Proof. If 9(t) is bounded, the number of integer values of k satisfying the 
second equation of or (*)2 for some (t, s) is necessarily finite. Then the 
result follows from Corollary [6j □ 

The converse of Theorem [7] is not necessarily true (see Example 2). So, we still 
need a characterization for the existence of infinitely many self-intersections. 

Lemma 8 The polynomial Kes s (a(t, s), (3(t, s, k)) cannot be constant. More 
precisely, it has positive degree in k. 

Proof. By Theorem [31 Res s (a(t, s), f3(t, s, k)) cannot be identically 0. More- 
over, by writing explicitly the system S = {a(t, s) = 0, (3(t, s, k) = 0}, 

A{t)B{s) - A(s)B(t) = 

< 

C(t)D(s) - C(s)D(t) - 2k7rD(t)D(s) = 

we observe that the points (£',£', 0), t! G R, are solutions of S. So, for any 
t' G R, Kes s (a(t, s), /3(t, s, k)) vanishes at (£', 0) and therefore it cannot be 
constant. For the same reason, Res s (a(t, s), /3(t, s, k)) cannot be an univariate 
polynomial in t. In fact, k is a divisor of Res s (a(t, s), /3(t,s,k)). □ 

Lemma 9 The polynomial Kes s (^(t, s), v(t, s, k)) cannot be constant. More- 
over, if(*2) has solutions for infinitely many values ofk G Z, then Res s (/i(t, s), u(t, s, k)) 
has positive degree in k. 

Proof. By Theorem [5l Res s (/x(t, s), u(t, s, k)) cannot be identically 0. More- 
over, by writing explicitly the system S' = s) = 0, u(t, s, k) = 0}, 

A(t)B(s) + A(s)B(t) = 

< 

C(t)D(s) - C(s)D(t) - (2k + l)-KD(t)D(s) = 

we observe that if A(t') = with D(t') ^ then there exists k' G K such that 
the point (t', t' , —1/2) is a solution of S. If D(t') = 0, then the points (t 1 , t', k) 
are solutions of S for any k. As a consequence, the polynomial Res s (fi(t, s), u(t, s, k)) 
cannot be a constant. 
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Now assume that (*) 2 has infinitely many solutions. That implies that there 
are infinitely (f , s', k') solutions of S with D{t') ^ 0, D(s') ^ 0, B{t') ^ and 
B(s') ^ . Due to Theorem O and the linearity of k in z/(t, s, k),the polynomial 
Res s (/i(t, s), u(t, s, k)) has positive degree in both t and k. □ 

Next we denote by £i(t,k) (resp. C,2{t,k)) the result of taking out from the 
square-free part of Res s (a(t, s), f3(t, s, k)) (resp. Res s (/i(t, s), v(t, s, k))) the 
univariate factors in t. 

Theorem 10 C has infinitely many self-intersections if and only if either 
£i(t,k) = or C,2(t,k) = is an algebraic curve (in the {t, k}-plane), non- 
bounded in k. 

Proof. If C has infinitely many self-intersections, by Corollary [6] one of the 
systems (*)i, (*)2 nas solutions for infinitely many values k £ Z. By reasoning 
over (*) 2 , one may see that each value t = t or s = s gives rise to 

just finitely many solutions of the system (otherwise, we would reach the 
conclusion that gcd(C, D) ^ 1). Consequently, if some system {*)\, (*)2 has 
solutions for infinitely many values k £ Z, S = {a(t, s) = 0,(3(t,s,k) = 0} 
or S' = {/J,(t, s) = 0, v(t, s, k) = 0} has infinitely many solutions too. Since 
resultants are combinations of the polynomials which define the systems S 
and S' and they are non-zero by Lemma [8] and Lemma [U either k) = 
or £ 2 (£j k) = is non-bounded in k. 

Conversely, assume that £i(t,k) is non-bounded in k (we would argue in a 
similar way with £ 2 (^ k) = 0). Since it has by definition no univariate factors 
depending on t, there are at most finitely many points (f, k') with k') = 
where the leading coefficient of a(t,s) with respect to s, D(t) and B(t) 
vanish. By the Extension Theorem of resultants, any other point of k) — 
corresponds to a solution of (*)i. Since £i(t, k) = is non-bounded in k, 
observe that it contains infinitely many points with k £ Z and so with t ^ s. 
Therefore C has infinitely many self-intersections. □ 

The above results are illustrated in the following examples. 

Example 1. Let C be parametrized in polar coordinates by ip(t) = (t,t). The 
plotting of this curve for t £ [— 5tt, 5tt] is shown in Figure |2J One may see that 
9{t) = t is not bounded; so, there might be infinitely many self- intersections. 
In this case we get k) = 2ktr and ^(t, k) = 2t — 2k7i — it. Finally, it is 
easy to see that k) = is bounded in k, but £ 2 (t, k) = is not. So, from 
Theorem [10] we conclude that C has infinitely many self- intersections. In fact, 
one may see that all these self-intersections lay on the y-axis. 
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Fig. 2. (p{t) = (t,t): Of) non-bounded, infinitely many self-intersections 



t 4 



Example 2. Let C be parametrized in polar coordinates by if if) = \ t, ^ 

T 3 3" 

Again 9 if) is not bounded. However, the plotting of this curve for t G — , - , 

shown in Figure El suggests that the curve has not infinitely many self- intersections 
(in fact, it has no self-intersections at all). Let us check this by using Theorem 
ITU1 We get £i(t, A;) = 2kix and £ 2 (t, = (2k + Both curves are clearly 
bounded in k; therefore, there are at most finitely many self-intersections. 




( t4 \ 

Fig. 3. ip if) = I t, ^ — - 1 : 0{t) non-bounded, finitely many self-intersections 

ft t 2 \ 

Example 3. Let C be parametrized by if if) = \ — — - , ^ 2 — - j . Since 9(t) = 
t 2 

is bounded, from Theorem [7] it follows that there are at most finitely 



t 2 + l 

many self-intersections. Furthermore, £ft, k) = 2kn{2kirt 2 + 2kir — t 2 + 1) and 
£ 2 (t, k) = ir(2k + l)(2/cvrt 2 + 2kix - t 2 + nt 2 + 1 + tt), both bounded in k. So, 
by Theorem [10] we derive the same conclusion. The plotting of the curve for 
t G (— oo, oo) is shown in Figure HI 
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3 7t 

2 



t t 2 



Fig. 4. <p(t) = ^-j : 9{t) bounded. 



£.5 Limit Circles, Limit Points and Spiral Branches 



Figure |5] shows the plotting of the curves defined by <pi(t) = 

for t E (0,6vr) (left), and tp 2 (t) = (t,l/t) for t E (0,tt/4) (right). In the 
first case, one sees that the curve winds infinitely around the circle r = 1, 
coming closer and closer to it. In the second case, the curve winds infinitely 
around the origin of coordinates, somehow "converging" to it. Finally, Figure 
[3] shows a curve that winds infinitely around the origin, but getting further 
and further from it. We will refer to the first situation by saying that here the 
curve exhibits a limit circle (r = 1 in this example). In the second case, we will 
say that the origin is a limit point; finally, we will refer to the third situation 
by saying that the curve presents a spiral branch. Along this subsection, we 
address these phenomena from a theoretical point of view and relate them to 
the appearance of infinitely many self-intersections. 



t^\ 

U 2 + i't 2 + i; 




Fig. 5. Example of Limit Circle (left) and Limit Point (right) 
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Definition 11 Let to G E U {±00}. We say that C exhibits, for t = to: 

(1) A limit circle if limt^ to r(t) = r G R, r 7^ 0, and limt_> to 9(t) = ±00. 

v4 limit point a£ £/ie origin if lim t ^ ) . to r{t) = 0, and limt^ to 9(t) = ±00. 
(3) A spiral branch^ if lim t ^t r(t) = ±00, and limt^ to 9(t) = ±00 

In eac/i case, we will say that to generates the limit circle, the limit point or 
the spiral branch. 

Notice that limit points are degenerated cases of limit circles, namely when 
ro = 0. One might wonder if limit circles can be centered at a point different 
from the origin, or if there can be limit points other than the origin. The 
answer, negative in both cases, is given by the following theorem. 

Theorem 12 A curve rational in polar coordinates cannot have a limit circle 
centered at a point different from the origin, or a limit point at a point different 
from the origin. 

Proof. If C had a limit circle centered at P 7^ (0, 0), or a limit point Q 7^ (0, 0), 
we would have infinitely many local maxima and minima of r(t). However, this 
cannot happen because r(t) is a rational function, and therefore the number 
of t- values fulfilling r'(t) = is finite. □ 

Remark 1 The same argument of Theorem[TM proves that a curve rational in 
polar coordinates cannot have any other "attractor" different from the origin, 
or a circle centered at the origin. 

Note that if C has infinitely many self-intersections, by Theorem [7]the function 
9{t) is not bounded and so there must be limit points, limit circles or spiral 
branches. More concretely, we have the following result. 

Proposition 13 Let I be a subset of E (not necessarily an interval) with 
infinitely many t-values generating self-intersections of C. Then the closure of 
I contains some t-value giving rise to either a limit point, or a limit circle, or 
a spiral branch of C. 

The converse of Proposition [13] is not true; for instance, in Example 2 we 
have a curve with a spiral branch for to = 00 without self-intersections. Now 
we want to characterize the situation when a limit point, a limit circle or a 
spiral branch have infinitely many close self-intersections. First we need the 
following definition. 

Definition 14 We say that a limit circle, a limit point or a spiral branch 
generated by t G E U ±00 has infinitely many close self-intersections, if there 



12 



exists some real interval I verifying: 

(1) t E 1; 

(2) I does not contain any other t-value generating a limit circle, limit point 
or spiral; 

(3) infinitely many t 6 / generate self-intersections of C. 

Definition [T^l is illustrated in Figure |6] the thin lines correspond to the branch 
generated by the interval / appearing in Definition [TH 




Fig. 6. Infinitely many close self- intersections 

Theorem 15 Let £ e R U ±oo generating a limit circle, a limit point or a 
spiral branch. Then 

(1) If to G R ; then to has infinitely many close self-intersections if and only 
if t = to is a horizontal asymptote of some of the curves k) = or 

6(t, k) = o. 

(2) If to = ±oo, then it has infinitely many close self-intersections if and 
only if some of the curves £i(t, k) = or £ 2 (t, k) — exhibits an infinite 
branch as t — > ±oo which is not a horizontal asymptote (i.e. there exists 
a sequence of real points (t n , k n ) of the curve with t n — > ±oo and k n — > 
±oo). 

Proof. We prove (1); the proof of (2) is similar. Assume that there exists 
an interval / C R satisfying Definition [TH Then the set of points of either 
£i(t, A;) = or ^{t.k) = with t 6 I must be non-bounded in k. This can 
only happen if there exists t a E I such that t — t a is an asymptote of either 
k) = or £ 2 (£, k) = 0. However, if t — t a is an asymptote then for any 
interval I a containing t a , and not containing any other t-value where 8(t) is 
infinite, we can find a non-bounded portion of either k) = or ^ 2 (^ k) = 0, 
therefore giving rise to infinitely many self-intersections of C, with t- values in 
I a . So, from Proposition [13] we have that t = t a . Using this same argument 
we can prove the converse statement, i.e. if t = to is an asymptote then any 
interval containing it has the desired properties. □ 
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The detection of asymptotes of an algebraic curve is addressed for example in 



Corollary 16 If t — t has infinitely many close self-intersections, then any 
interval containing t generates infinitely many self-intersections of C. 

The above results are illustrated in the following examples. 

Example 1 (cont.): Recall that £i(t, k) = 2kir and £ 2 (£, k) = 2t — 2kix — ix. 
The second one has an infinite branch as t — > ±00; so, from Theorem [TBI we 
deduce that every t-interval of the form (—00, a) or (b, 00) contains infinitely 
many t- values generating self-intersections of C. 



curve has a spiral branch for t = and a limit point for t — > ±00. Observe 
Figure [7| where the curve is plotted for different values of t. A direct compu- 
tation yields k) = k{kixt + 1) and 

&(t,k) = 4t 6 -4vr (2 k+ l)t A -4t 3 + tt 2 {2k + if t 2 + 2tt {2k + l)t + 2. 



The curve ^{t, k) = is empty over the reals. However a factor of k) = 
corresponds to a hyperbola whose asymptotes are t = and k = 0. Hence, 
from statement (1) of Theorem [T5l we deduce that the spiral branch generated 
by £o = has infinitely many self-intersections. 



Example 4. Consider the curve defined 





t 6(1,5) 



t 6 (-5,-1) 



t 6 (-5,-0.15) U (0.15,5) 
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2.4 Asymptotes 



It is classical (see for example [5]) that asymptotes in polar form correspond 
to values t = t G K U {±00} such that: 

• \im t ^ to r(t) = ±00. 

• lim t _> to 0(t) = «6R. 

• Hm t ^ to r(t) • (0(t) - a) E 1. 

If the above conditions hold, then the asymptote is the line parallel to the line 
with slope a, at distance 5 = \im t ^t r(t) ■ (9(t) — a) of the origin. 



lim t _> ±00 r(t) = ±00, lim^ioo 9(t) = 1 and \im t _, ±00 r(t) ■ (9(t) - 1) = 0. So, 
the line passing through the origin and with slope 1, i.e. y = x, is an asymptote 
of the curve for t — > ±00 (see Figure [8]). 



3 Visualizing Curves that are Rational in Polar Coordinates 

The goal of this section is to introduce an algorithm, polares, for visualizing 
the "interesting" part of C, providing as well relevant information about the 
curve. We have implemented it in Maple 15 and tested it over several examples, 
some of which are presented in this section. 

Next we describe the algorithm in the following different cases: 

(1) r(t) and 6(t) are both bounded 

(2) 6{t) is bounded and r{t) is not. 

(3) r(t) is bounded and 6{t) is not. 



Example 5. Consider the curve 




Then we have that 




Fig. 8. An asymptote 
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(4) r(t) and 9(t) are both unbounded. 

Obviously, the first step of the algorithm is to detect the case we are in; then 
the algorithm proceeds accordingly. 



3.1 Algorithm polares when r(t) and 6(t) bounded 

Input: A proper polar parametrization (p(t) = (r(t),8(t)). 
Output : 

(1) Information about the existence of P^. 

(2) Information about the self-intersections. 

(3) Plot of ip(t) for t in M using the Maple command polarplot. 
Examples: 

• v>i(t) = (tw>tw) 

>polares([ T ^, T ^]); 
r and theta both bounded 

Real point at the infinity such that (r, theta) =[0, 1] and the point 
is [0, 0] 

The point at infinity (0,0) is reached 1 times in R, so self -intersection 
at the origen 



\+t 2 ' i+t 2 
t t 2 +u 



3 Jt X^'^ 
4 ^ 




ff><: 






\0.1 0.2 q'.3 W 

31/ 

4 


3 

plott=- 


00 00 


Fig. 9. <p x {t) 



>polares([ T ^,\^]); 
r and theta both bounded 

Real point at the infinity such that (r, theta) =[0, 1] and the point 
is [0, 0] 



1(3 



The point at infinity (0,0) is reached 1 times in R, so self -intersection 
at the origen 

System (1) gives self -intersections for k in [[ -2,2]], k<>0 
System (2) gives self -intersections for k in [[ -2,1]] 




2 

plot t=- oo ..oo 



Fig. 10. <f>2(t) 

3.2 Algorithm polares when 9 (t) is bounded and r{t) not 

Input: A proper polar parametrization ip(t) = (r(t),0(t)). 
Output : 

(1) Information about the existence of Poo. 

(2) Information about the self-intersections. 

(3) Information about the existence of asymptotes. 

(4) Plot of (p(t) using the Maple command polarplot: 

(a) If there are not values of t generating asymptotes, then we plot the 
curve for t G (— oo, oo); 

(b) Otherwise, let Ti be the set of values of t generating asymptotes, let 
T 2 be the set of real values of t such that r(t) = and, if \9(t)\ < 2 n 
for t G K, let T 3 be the real values of t generating the maximum and 
minimum of 0(t). 

Let P : = Ti U T 2 U T 3 = {t u . . . , t m } with h < . . . < t m . 
Now, 

(i) If {ti,t m } ^ {oo, — oo}, then G M 2 and its plot may be of 
interest. So we add t := — oo and t m+ i := oo to the set P. 
For every ti G P, 1 < % < m, we plot the curve tp(t) for t G 
(^pi, **±fs±!). For t in ( tj± ^-, U), the color of the plot will be 
red and for (ij, ) the color will be blue. 
As for ti and t m , we plot <p{t) in the ranges {2t x — ^y 2 -, ^f 1 ) 
and ( t "' 2 tro - 1 , 2t m - 

As for t and t m+ i, we plot ip(t) in the ranges (— oo, 2t\ — ^f 2 -) 
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and (2t m - tm+ * m+1 ,oo). 
(ii) If {ti,t m } = {oo, — oo}, then for every ti £ P, 2 < i < m — 1, 
we plot the curve tp(t) in the range ( t ' + *'~ 1 , tl+ *' +1 ). For i in 
( * + 2 i ~ 1 , U), the color of the plot will be red and for (U, ti+t ^ +1 ) 
the color will be blue. 

As for t 2 and t m _i, we plot ip{t) for the ranges (t 2 — 10, ^rp 3 -) 
and (*s=i=*2=a,t m _ 1 + l0). 

Finally, as for {oo, — oo}, we plot y>(£) for the ranges (£2 — 20, — 

10) and (t m _i + 10,t m _i + 20). 
Let us point out that we border the values of t generating the asymtotes 
for plotting. 

Examples: 

• <Pa(t) = ( t 2_i lt+30 , t -j0) 

>polares ( [t 2 /(t 2 - 11 1 + 30), (t 2 + 78)/(t 2 + 1)] ) ; 
r unbounded and theta bounded 

Real point at the infinity such that (r, theta) =[1, 1] and the point 
is [cos(l), sin(l)] 

Point at infinity is not reached with k=0 
Point at infinity is not reached with k<>0 

System (1) gives self-intersections for k in [[ -2,2]], k<>0 
The values of t generating asymptotes are {5., 6.} 
Values of t considered in the plot { — oo, 0., 5., 6.00} 



<Pi(t) = (t, 



1+t 2 ' 



>polares ( [t, (t 2 + 14)/(t 2 + 1)] ) ; 

r unbounded and theta bounded 

There is no point at infinity 

Values of t generating asymptotes [00, —00] 

Values of t considered in the plot {0.,oo,— 00} 



3.3 Algorithm polares when r(t) is bounded and 6{t) not 



Input: A proper polar parametrization ip(t) = (r(t),0(t)). 
Output : 

(1) Information about the existence of P^. 

(2) Information about the existence of limit circles. 

(3) Information about the existence of limit points. 

(4) Information about the self-intersections. 
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plot around the value of t equal to -infinity 



plot around the value of t equal to 



plot around the value of t equal to 5 




plot around the value of t equal to 6 plot around the value of t equal to infinity 

Fig. 11. <ps(t) 




plot around the value of t equal to - °o plot around the value of t equal to 0. plot around the value of t equal to « 

Fig. 12. <p A (t) 

(5) Plot of ip(t) for t in K using the Maple command polarplot. 

Let Ti be the set of values of t generating limit circles, let T 2 be 
the set of values of t generating limit points, let T3 be the set of 
real values of t such that r(t) = and let T 4 be the real values of t 
generating the maximum of 

Let P := 71 U T 2 U T 3 U T 4 = {h, . . . , t m } with t x < . . . < t m . Then 
we proceed similar to the case (BJ of Section 13.21 Now we border the 
values of t generating limit circles for plotting. 

Example: 

. ^ 5 (t) = (t 2 /(* 2 + i),t7(* 2 + i)) 
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>polares ( [t 2 /(t 2 + 1), t 3 /(t 2 + 1)] ) ; 

r bounded and theta unbounded 

There is no point at infinity 

Values of t generating limit circles [—00,00] 

There are infinitely many self-intersections 

t=infinity has infinitely many close self-intersections 

There are no limit points 

Values of t considered in the plot {O.,oo,— 00} 




plot around the value of t equal to - » plot around the value of t equal to 0. plot around the value of t equal to °° 

Fig. 13. <p 5 (t) 



3.4 Algorithm polares when both r(t) and 9{t) are unbounded 

Input: A proper polar parametrization ip(t) = (r(t),6(t)). 
Output : 

(1) Information about the existence of P^. 

(2) Information about the existence of limit circles. 

(3) Information about the existence of limit points. 

(4) Information about the existence of spiral branches. 

(5) Information about the existence of asymptotes. 

(6) Information about the self-intersections. 

(7) Plot of tp(t) for t in K using the Maple command polarplot. 

Let Xi be the set of values of t generating limit circles, let T 2 be the 
set of values of t generating limit points, let T 3 be the set of values of 
t generating spiral branches, let T 4 be the set of values of t generating 
asymptotes and let T 5 be real the values of t such that r(t) = 
Let T 1 UT 2 UT 3 ur 4 UT 5 = {t u . . . , t m } with t x < . . . < t m . 
Then we proceed similar to the case (TjQ) of Section I3.2L Here we 
border the values of t generating asymptotes, limit circles and spiral 
branches for plotting. 
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Example: 



• <p 6 (t) = (t,(t 3 + l)/(t 2 -3t + 2)) 

>polares ( [t, (t 3 + 1) /(t 2 - 3 1 + 2)] ) ; 

r and theta both unbounded 

There is no point at infinity 

Values of t generating limit circles [1., 2.] 

There are no limit points 

Values of t generating spiral branches [—00,00] 
There are not values of t generating asymptotes 
There are infinitely many self-intersections 
t=l has infinitely many close self-intersections 
t=2 has infinitely many close self-intersections 
t=infinity has infinitely many close self-intersections 
Values of t considered in the plot {— 00, 0., 1., 2., 00} 




plot around the value of t equal to 2. plot around the value of t equal to °o 



Fig. 14. <p 6 (t) 
4 Conclusions/Further Work 

In this paper we have presented several theoretical results and an algorithm 
for properly plotting curves parametrized by rational functions in polar form. 
Our results allow to algorithmically identify phenomena which are typical 
of these curves, like the existence of infinitely many self-intersections, spiral 
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branches, limit points or limit circles. Furthermore, the algorithm has been 
implemented in Maple 15, and provides good results. Natural extensions of 
the study here are space curves which are rational in spherical or cylindrical 
coordinates, curves which are algebraic, although non-necessarily rational, in 
polar coordinates (i.e. fulfilling h(r,6) = 0, with h algebraic), and similar 
phenomena for the case of surfaces. It would be also interesting to analyze the 
curves defined by (implicit) expressions of the type f(r, sin(0), cos(0)) = 0, 
where / is algebraic, since this class contains, and in fact extends, the class 
of algebraic curves; also, it includes the important subclass of curves defined 
by equations r n = g(9), with g(9) a rational function, which often appear in 
Geometry and Physics. Some of these questions will be explored in the future. 
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